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Abstract: Leptoquarks in the mass range of 500-1000 GeV can be resonantly produced 
in significant numbers by PeV neutrino interacting with nuclei at IceCube. We compute 
the event rates of leptoquark production and decay events and use the 3-year IceCube 
data for PeV energy events to find the allowed range of the leptoquarks mass and coupling 
parameter space. We use a low-scale quark lepton unification model based on the SU (4)c(8> 
SU{2)l ® U{l)ii gauge group where leptoquark couplings which give rise to proton decay 
are forbidden by the symmetry. We constrain the parameters of this model and point out 
signals of leptoquarks in this model which may be seen in PeV energy IceCube events in 
the future. 
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1. Introduction 

Leptoquarks (LQs) are SU{3) coloured particles which also have non-zero lepton as well 
as baryon numbers. These arise in grand unified theories like SO(IO), SU(5) or Pati-Salam 
SU(4) and are expected to have masses in the GUT scale in order to not cause rapid proton 
decay. Models with low scale leptoquarks have been considered [1-4] which do not lead to 
rapid proton decays. The leptoquark model [4] is based on the Glq = SU{4:)c 'S' SU{2) 1 0 
U (1)_R quark-lepton unification group which breaks to the standard model at some low scale 
~ TeV. The SM quarks and leptons in this model can be unified in the same multiplet 
(QlJl) ~ (4,2,0), {uR, Ur) ~ (4,1,-1/2) and {dR,eR) ~ (4,1,1/2). The scalars of this 
model give rise to this symmetry breaking and also to provide Yukawa couplings to generate 
fermion masses. They do not couple to the type of fermion-bilinears [5](Qx,Ti, Ql'Ql etc.) 
which give rise to proton decay. In addition to these scalar leptoquarks there are vector 
leptoquarks from the gauge bosons of the Gqr gauge group but these are constrained to 
be heavy (Mx > 10^ TeV) to evade bounds from rare meson decays [6,7]. 

Leptoquark can be produced, singly or doubly, at colliders (hadron, e'’'e“ or e^p) and 
signals (like the Ujj from the decay of a leptoquark pair) and other important properties 
and constraints can be studied [8-19]. A search for pair-production of first and second 
generation scalar LQs has been performed with 19.6 fb“^ of data by CMS [20,21]. According 
to this the first generation scalar LQs with masses less than 1005 (845) GeV are excluded 
for (5 = 1(0.5), where /3 represents the branching fraction of an LQ to a charged lepton 
and a quark. For second generation scalar LQs the exclusion limits are 1070 (785) GeV 
for /3 = 1(0.5). Similarly, scalar LQs of third generation with masses below 740 GeV are 
excluded at 95% C.L. assuming a 100% branching fraction for the LQ decay to a tau lepton 
and a bottom quark [22]. A recent 2.6c7 excess in a iijj$T events, reported by CMS [23], 
has been explained by various models of LQ in the mass range 550-650 GeV [24,25]. 

It has been pointed out [26] that a natural arena for the production of leptoquarks is 
the neutrino-nucleon interactions at IceCube [27]. Recently the leptoquark interaction of 
the form S^{£l.Ql + u^r^) of a scalar leptoquark with the third generation leptons and 
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first generation quarks has been proposed as an explanation of the 2 year IceCube data [28]. 
A very recent account of a minimal leptoquark scenario in the light of CMS and IceCube 
observations are presented in [29]. Restriction to the third generation leptons explains the 
paucity of muon tracks at neutrino energies 0.1- 2 PeV and the cross section is enhanced 
by the resonant production of ~ 600 GeV leptoquarks at neutrino energies ~ PeV. 

In this paper we discuss the production and decay of the scalar leptoquarks of the low 
scale quark-lepton unification model based on SU{4)c ® SU{2 )l ® U{1)r [4] at IceCube. 
The leptoquarks in this model decay into hadronic showers (plus neutrinos). The analysis 
of the 988 days data gives a neutrino flavor ratio of 1:1:1 consistent with all 37 events 
being from standard neutrino CC and NC events [30]. Thus to put an upper bound on 
the parameters of the leptoquark model we make the conservative assumption that no 
event among the 35 with less than PeV energies are of leptoquark origin. The resonant 
production of leptoquarks with masses in the 500-1000 GeV, however, can give significant 
number of hadronic shower events at IceCube at PeV range. 

2. Quark-Lepton Unification Model 

We have used the low-scale quark-lepton unification model described in [ 4 ] which predicts 
the existence of both vector and scalar leptoquarks . This leptoquark model is based on 
the Gql = SU ( 4 )c ( 8 > SU {2)l®U{1)r quark-lepton unification group which spontaneously 
breaks to the standard model Gsm = SU{3) 0 SU{2)l (g) U{1)y at some low scale. The 
SM quarks and leptons can be unified in the same multiplet {Ql-,^l) ~ (4, 2,0), {ur, vr) ~ 
( 4 , 1, —1/2) and {dR, cr) ~ ( 4 , 1, —1/2). We concentrate on the scalar LQs of this theory, 
namely, 4*3 (3,2, -1/6)gj^ and 4>4 (3,2,7/6)gj^ which are present in the scalar 4> of 

the theory. The masses of the components of <^3 and 4>4 are determined by the parameters 
of the scalar potential. The relevant part of the potential is given by [4], 

V = m|Tr[4>t4>] -F A2i7^'i7Tr[4>t4>] -k A3X^yTr[4>t4>] -k 

+ A9Tr[$t^,]2 + Aio(Tr[$t$]) 2 ^ ( 2 . 1 ) 

where under Grq, ~ (15, 2 , 1 / 2 ), H ~ ( 1 , 2 , 1 / 2 ) and y ~ (4, 1 , 1 / 2 ). Also, 

where is the SU{A)c generator and H 2 is another Higgs doublet. Moreover the vacuum 
expectation values (vev) of the scalar fields are given by {x^) = v^ly/2, (H^) = ui/\/2 and 
{H 2 ) = U2/\/2. With these vevs the mass terms for the fields ‘I> 3 ^ 4^8 can be written as, 

(^m| + ^ (^(A 2 + A5)u? + + ^4^4 + ^^4*8) + (2.3) 

From the mass terms it is evident that 4>8 and 4>4 are degenerate in mass but <^3 mass 
can be different from the former two, owing to the Ae term. In this context it is worth 
mentioning that the LHC searches at ^/s = 8 TeV with 19.7 fb“^ of integrated luminosity. 
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through a dijet resonance above 1.2 TeV have excluded the mass range of [1.3-2.5] TeV of 
the color-octet mass [31]; also according to a more recent search the bound on this mass 
exists down to 500 GeV [32]. Thus the octet <1>8 and the leptoquark <^4 can be well within 
this bound but the leptoquark <^3 can be of smaller mass by taking small negative Ag¬ 
in this work we take them to be the free parameters which we constrain from the 
IceCube data. The relevant interactions of the leptoquarks are, 

>Clq = A 2 QL‘h 3 i^*^ + \ 2 ^l^avP + -|- h.c. (2.4) 

The symmetries of LQs $3 and $4 do not allow interaction terms with type of fermion 
bilinears {Ql ■ ^l,Ql ■ Ql) that lead to proton decay via these LQ exchanges at tree level 
at dimension four. However, there exists dimension five operators that may lead to proton 
decay. But these dimension five operators can be avoided by the imposition of appropriate 
discrete symmetries [5]. The constraints on these types of LQs coming from charged lepton 
sector, e.g., via the processes jjL ey, conversion of /r to e and electric dipole moment of 
the electron, have been discussed in [5]. 

The coupling A 4 < 0.01 and this constraint comes from the kaon decay —)• 

to be in agreement with fermion masses [7]. The coupling A 2 is not constrained and we 
put bounds on this coupling from IceCube events. Constraints on leptoquark models from 
other rare decays has been considered in [33,34]. 


Leptoquark 

Electromagnetic Charge 

(Q) 

Neutrino coupling 

4*3 

+ 2/3 

X2i'u 

-1/3 

X2vd 

4?4 

+5/3 

- 

+2/3 

X2i'u 


Table 1: Leptoquarks with appropriate electromagnetic charges and couplings. 

In Table 1 we present the electric charges of the components of the LQs and their 
couplings with the neutrinos. Note that the state of $4 with electric charge (-1-5/3) will 
have no coupling with neutrinos. All the processes relevant for the production of LQ from 
neutrino-nucleon interaction and its subsequent decay to the r and hadronic showers are 
(here the superscripts in < 1 > 3^4 represents the corresponding electric charges.), 

vu —>• —>■ T'u, i>d —)• $3 —)• vd, i>u —>• —>■ T'u , (2.5) 

these processes are proportional to the coupling A 2 on which the constraints are mild. 

On the other hand the model (2.4) allows the following processes which could have 
contributed to the IceCube signal, 

uu —> <I> 3 ^^ — i'^d, vd —> 4)3 > l~u, uu —> <I> 4 ^^ —> £~^d ( 2 . 6 ) 
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but these processes are proportional to the coupling A 4 which is already constrained to be 
A 4 < 0.01 from kaon decay —)• [7]. We study LQ production and decay processes 

where only the A 2 coupling arises. Due to the smallness of A 4 the LQ decays to + j, 
which is the signal in [28], is not significant in this model. We considered A 2 to be 0.5 and 
1.0 and we keep A 4 to a fixed value 0.005 in our subsequent study. 

The neutrino-nucleon cross section for the production and subsequent decay of LQ of 
mass M^), in the narrow width approximation, is given by [35,36], 

cf'LQivN 9X) = BR(^ ^ 9q)xqNix-n‘^), (2.7) 

where qN{x] /i^) is the parton distribution function (PDF) of the parton q in the nucleuon 
(with mass mjv) and x = /{2mNEy) is the parton fractional momentum. The renor¬ 
malization scale fi is set at M$. We have used the nuclear PDFs, CTEQ611 [37] at LO to 
get the qN{x\ jjp). 



(a) 



(c) 




Figure 1: The distribution of neutrino-nucleon cross section for various cases. The values 
of the coupling A 2 that have been used are shown in each plot. 


In Fig. 1 we show the distribution of neutrino-nucleon cross section, —>• i^X). 

The upper panel (Figs, la and lb) shows the case where only the leptoquark $3 ' is 

9/0 

involved. In the lower panel (Figs. Ic and Id) we show the distribution where either ‘kg' 

2/3 

or $ 4 ' comes in the intermediate state. 
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3. Constraints from IceCube Events 


The three-year (2010-2012) data set, with a livetime of 988 days, reveals 37 neutrino events 
in the energy range 0.3-2 PeV. Among these, 9 are track events and 28 are cascade events, 
with the flavor ratio 1:1:1 [30] expected from pion/muon decay neutrinos oscillating 
over cosmological distances [38,39]. The astrophysical neutrino flux (averaged over zenith 
angle) follows a power law [27], 


/ f? \ ~2.3 

$(EJ = 4.74 X 10"^ GeVcm-^s-^sr-^ (3.1) 

\lGeV / 

The highly energetic incoming neutrinos can produce leptoquarks after the collision 
with quarks in the nucleons at IceCube. The subsequent decay of the leptoquarks will be 
registered as events at IceCube. The expected number of events at IceCube is given by. 


A/lq = j dE^[ai,Q{E^).Bv]^{Eu), (3.2) 

where ut is the effective number of target nucleons in IceCube and is ~ 6.0 X 10^^ as the 
effective volume is roughly 1 km^; the time of exposure T = 988 days. (Tlq(£'j/) has been 
defined in Eq. 2.7; the Br represents the branching ratio of the leptoquark d* decaying to 
a neutrino and a quark. 

The event rate distributions, are shown in Fig. 2. It is worth-mentioning that the 
rapidly falling nature of these distributions are dictated by the spectrum of neutrino flux 
HE,). 

To test the allowed parameter space of the A 2 coupling and LQ mass for extra LQ events 
in the IceCube observations we do a y-square fit with {X 2 , M 1 / 3 ) as free parameters^. 
We define. 


X^{H, M -1/3) 

^3 


E 

energy bins 


(A/" ^xp 


+ A/lq))^ 
A2 


(3.3) 


where the error A is the experimental error [27] in the measured events in each energy bin 
and the uncertainty in the model prediction added in quadrature. The events expected 
from the standard model CC and NC interactions from atmospheric neutrinos plus the 
extra-terrestrial neutrinos with flux spectrum (3.1) have been computed in [27]. We add 
to the SM events the events expected from LQ processes and minimize the summing 
over energy bins from E, = 100 GeV — 10 PeV. The minimum is achieved for A 2 = 0 
which means that LQ events spoil the overall fit from SM events. When the LQ mass is 
M i/3 < 500 GeV then it gives an excess contribution over observations in the E, < 1 

®3 

PeV energy bins where the fit from standard model CC and NC events from atmospheric 
neutrinos and extra-terrestrial neutrinos with spectrum (3.1) fits observed IceCube events 

^Similar studies can be made for H/t also. 
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(a) (b) 




(c) (d) 

Figure 2: The event rate distribution dN/dE^, for various cases. The coupling A 2 is shown 
in each plot. 


well. When the LQ mass is M^- 1/3 > 650 GeV (and the corresponding A 2 > 1, in order to 
maintain the cross section) there are excess events in the ~ 6 PeV Glashow-resonance 
region where there is already a problem in that, even the SM prediction is not observed. 
In Fig. 3 we show the contribution of LQ with mass M 1/3 = 500 GeV and M 1/3 = 650 
GeV respectively to the IceCube events. We see that in the mass range M^- 1/3 = 500 — 650 
GeV with coupling A 2 = 1 there is an improvement in the fit for the events in the 1 — 2 
PeV bins compared to the standard model. 

As the total number spectrum of events has a worse fit with the inclusion of lepto- 
quark events, we can rule out some parameter space of LQ models from IceCube events. 
In Fig. 4 we show the (A 2 ,M^_i/ 3 ) parameter space where the region below the curve is 
disallowed at 95% CL. 

4. Summary and Conclusion 

Leptoquarks in the mass range 500-1000 TeV can be produced in significant numbers by 
the PeV neutrinos at IceCube. Assuming that the 37 events seen in the IceCube data can 
be explained by t'e, t'/x, CC and NC events we put bounds on the masses and couplings 
of scalar leptoquarks of the low scale quark-lepton unification model based on the gauge 
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Figure 3: Event rate distribution for leptoquark with A 2 = 1 and M^- 1/3 = 500 GeV (left) 
650 GeV (right). 



-^2 


Figure 4: The (A 2 , M^- 1 / 3 ) parameter space where the region below the curve is disallowed 
at 95% CL. The Xmin = 3.66 for A 2 = 0 which means that when all the energy bins are 
considered in the the SM gives the best fit to the observed data. 


group SU{A)c ® SU{ 2 )l ® U{1)r [4]. Of the possible signals in this model the vq —)• i)q 
process which will give rise to a hadronic shower that can be observed at IceCube. The 
9q —)• Iq process allowed in the model which would produce lepton and hadronic showers 
is constrained to be small from collider bounds on the corresponding coupling. We put 
bounds on the scalar leptoquark masses which connect the first generation quarks with 
the third generation leptons to be above the 500-1000 GeV if the coupling A 2 is in the 
0.1-1 range. Observations with longer exposures over target volumes can produce hadronic 
showers by resonant leptoquark decay which may be observed over the background. 

In the context of colliders, we would like to add that in this model the eejj or evjj 
signal, studied in [19,25,29], will take place only via coupling A 4 (see Eq. 2.4), but this 
coupling is restricted by rare meson decay processes, which constrain A 4 < 0.01 [7]. For 
these values of A 4 the leptoquark events at IceCube are negligible compared to those gen¬ 
erated by A 2 coupling. This model does not explain the eejj or evjj signals at LHC and 
is only constrained by IceCube data. Also due to the very nature of the couplings, in this 
model, one can have vvjj type of signal but that will be difficult to examine at colliders. 
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